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Abstract 

The linear-fractional Galton-Watson processes is a well known case when many characteristics of a 
branching process can be computed explicitly. In this paper we extend the two-parameter linear-fractional 
family to a much richer four-parameter family of reproduction laws. The corresponding Galton-Watson 
processes also allow for explicit calculations, now with possibility for infinite mean, or even infinite 
number of offspring. We study the properties of this special family of branching processes, and show, 
in particular, that in some explosive cases the time to explosion can be approximated by the Gumbel 
distribution. | 

1 Introduction 

Consider a Galton-Watson process (Z n ) n >o with Zq = 1 and the offspring number distribution 

Pk = P(Zi = k), k> 0. 

The properties of this branching process are studied in terms of the probability generating function 

f(s) =Po+Pis+ P 2 S 2 + ..., 

where it is usual to assume that /(1) = 1, however, in this paper we allow for /(1) < 1 so that a given 
particle may explode with probability p^ = 1 — /(1). The probability generating function f n (s) = E(s Zn ) 
of the size of the n-th generation is given by the n-fold iteration of f(s) 

/o(s) = s, f n {s) = /(/„-i(s)), n > 1, 

and therefore it is desirable to have a range of probability generating functions / whose iterations can be 
computed explicitly. 

The best known case of explicit calculations is the family of linear-fractional Galton-Watson processes 
with 

f(s)=Po + (l-Po), -rp— v-j [0,(1-p)" 1 ), 

1 - (1 -p)s 

representing the family of modified geometric distributions 

Pk = (1 -Po)(l ~p) k ~ X P, k> 1, 

fully characterized by just two parameters: po € [0,1) and p € (0,1]. In Section [2] for each 9 £ [—1,1] 
we introduce a family Qg of functions with explicit iterations containing the linear-fractional family as a 
particular case. In Section [3] we demonstrate that all f £ Qg are probability generating functions with 
/(1) < 1. A Galton-Watson processes with the reproduction law whose probability generating function 
belongs to Qg will be called a theta-branching process. 

The basic properties of the theta-branching processes are summarized in Section |4] where it is shown 
that this family is wide enough to include the cases of infinite variance, infinite mean, and even non-regular 
branching processes with explosive particles. 
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Recall that the basic classification of the Galton-Watson processes refers to the mean offspring number 
m = EZ\. Let q € [0,1] be the smallest non-negative root of the equation f(x) = x and denote by 

T 0 = inf{n : Z n = 0} 

the extinction time of the branching process. Then q = P(Tq < oo) gives the probability of ultimate 
extinction. For m < 1 and p\ < 1, the extinction probability is q = 1, while in the supercritical case m > 1, 
we have q < 1. 

If /(l) < 1, then the Galton-Watson process is a Markov chain with two absorption states {0} and {oo}. 
In this case the branching process either goes extinct at time T 0 or explodes at the time 

Ti = inf{n : Z n = oo}, 

with 

P(Ti < n) = 1 — /„(1), P(Ti < oo) = 1 - q, 

where the latter equality is due to /„( 1) — > q. In Section[5j using explicit formulas for f n (s) we compute the 
distribution of the absorption time 

T = To A Ti. 

Note that in the regular case, we have P(T\ = oo) = 1 and therefore, T = Tq. Observe also that the case 
/(1) < 1 has other, biologically more relevant interpretations. For example in the multitype setting, T\ can 
be viewed as the time of the first mutation event, see [9] . 

Also in Section [5] we consider a situation when the explosion of a single particle has a small probability, 
so that T\ takes large values in explosion scenarios. We show that in such a case the time to explosion can 
be asymptotically characterized with help of a Gumbel distribution. In Section [5] we study the Q-processes 
for the theta-branching processes extending the classical definition to the non-regular case. Our explicit 
calculations demonstrate that in the non-regular case the behavior of a branching process is more similar 
to that of the subcritical rather than supercritical regular case. Using these results on the Q-processes we 
derive the conditional limits of the theta-branching processes conditioned on non-absorption. 

A remarkable property of the linear-fractional Galton-Watson processes is that they can be embedded 
into the linear birth-death processes. In Section [7] we establish embeddability of theta-branching processes. 


2 Probability generating functions for theta-branching processes 


Using an alternative parametrization for the linear-fractional probability generating functions, we obtain 


1 


a 


1 — s 


+ c, 


se [0,1), 


(i) 


where 


P 1 ~P 

a = --, c = --. 

1 - Po 1 - Po 

This observation immediately implies that the n-fold iteration f n of the linear-fractional / is also linear- 
fractional 


1 


+ c(l -f Q ...-{- d r ‘ 
The key idea of this paper is to expand the family © by 


1 - f n (s) 1 - s 


{A — f(s )) s — a(A — s) 9 + c, s€ [0,i), 
with the help of two extra parameters (A, 9) which are invariant under iterations. 


( 2 ) 


Definition 1 Let 9 £ (—1,0) U (0,1]. We say that a probability generating function f belongs to the family 

Qe if 

f(s) = A — [o(A — s)~ e + cy 1 / 0 , 0 < s < A, 
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where one of the following three options holds 


(f) 

a > 1, 

c > 0, 

0£ (0,1], 

A = 1, 

(u) 

a £ (0,1), 

c=(l- a){ 1 - q)~ 9 , 

q£ [0,1), 

A = l, 

(Hi) 

a £ (0,1), 

c = (1 - a){A- q)~ 9 , 

q e [o,i], 

A>1. 

Definition |T] can be extended to the case 

9 = 0 by the following 

continuity argument 


A - [a(A - s)~ 9 + (1 - a)(A - q)- 9 }- 1 / 9 -► A - {A - g) 1 "^ - s ) a , 0^0. 
Definition 2 We say a probability generating function f belongs to 

• the family Go if for some a £ (0,1), 

f(s) = A-{A-q) 1 - a (A-s) a , 0<a<A, 

where either A = 1, q £ [0,1), or A > 1, q £ [0,1], 

• the family f € G-i if for some q £ [0,1] and a £ (0,1), 

f(s) = as + (1 — a)q , 0 < s < oo. 


Definition 3 A Galton- Watson process with the reproduction law whose probability generating function 
f £ Qe, 8 £ [—1,1], will be called a theta-branching process. 

It is straightforward to see, cf. Section [4j that each of the families Ge is invariant under iterations: if 
f £ Go, then /„ £ Ge for all n > 1. The fact, that the functional families in Definitions |T| and [2] are indeed 
consist of probability generating functions with /(1) < 1, is verified in Section [3] 

Parts of the Ge families were mentioned earlier in the literature as examples of probability generating 
functions with explicit iterations. Clearly, G i U G-i is the family of linear-fractional probability generating 
functions. Examples in [12; leads to the case A = 1 and 9 £ [0,1), which was later given among other 
examples in Ch. 1.8 of m- The case A = 1 and 9 £ (0,1) was later studied in IT'. A special pdf with 
8 = — 1 / 2 , 

/(a) = 1 - {ay/ 1 - s + 1 - a) 2 , a £ (0,1), 

can be found in [3] on page 112, as an example of non-regular Galton-Watson processes. 

Notice that there is a version of linear-fractional Galton-Watson processes with countably many types 
of particles, see [7]. It is an open problem to expand the theta-branching processes with 9 £ (—1,1) to the 
multitype setting. 


3 Monotonicity properties 


It is straightforward to see that each / € Go is a probability generating function with 
f'( s ) = (A — q) 1 ~ a a{A — s) a_1 , 

/ (n) (s) = {A - q) 1 ~ a a{ 1 - a)... (n - 1 - a){A - s) a ~ n , n > 2, 

and 


Po = A — {A - g) 1 -^, 

Vl = {A-q) 1 - a aA a -\ 
n — a — 1 

p n =p n _ i-— , n > 2. 

Therefore, {p n )n>i are monotonely decreasing with 

n 1 

p n = aA a (A - q) 1 ~ a A~ n (l - , 

fc =2 


so that p n ~ const •A n n 1 “asn->oo. 


n > 2, 
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Proposition 4 Let 6 £ (—1, 0)U (0,1) and f £ Ge- Then f is a probability generating function with /(1) < 1 
such that 


Po = A - (aA 9 + c) 1/9 , 
p^aia + cA 9 )- 1 - 1 / 9 , 


and for n > 2, 


aA~ n+1 


y ( cA " ) 

A- 1 V a + cA e ) 


Pn = 


Bi 


(a + cA e ) i ^n! ~'y a + cA6 ' 
where all B i n = B i n (6) are non-negative and, for n > 2, satisfy the recursion 


Bi^n — (rt 2 i6)Bi tn -i + (1 + ^ 0) LI^ — 1 5 72—17 ^ — lj..., n 1, 
with Bo n = B n ^ n = 0 for n > 1, and B\ 2 = 1 + 0. 

PROOF In terms of 

m ■■= = [a + c(A - s) 9 }- 1 / 9 , 0'(s) = c(A - s) e -V(s) 1+e , 

we have 


f{s) = acf>(s) 1+e , 

f"(s) = (1 + 0)ac(A - s) 6 ~ V(s) 1+2e , 

f"{s) = (1 + 0)(1 - d)ac(A - s) 9 ~ 2 (/){s) 1+29 + (1 + 0)(1 + 2 d)ac 2 (A - s) 29 ~ 2 (j)(s) 1+3e . 


and more generally, 

n —1 

f (n \s) = 5Z B ^ac\A - s) i0 -” + V(s) 1+ ( i+1)e , n > 2, 

i=1 

where Bi iTl are defined in the statement. To finish the proof it remains to apply the equality p n = f( n \0)/n\. 


In the linear-fractional case we have pk > Pk+i for all k > 1. The next extension of this monotonicity 
property was first established in El- 

Corollary 5 Let 6 £ (0,1) and f £ Gs with A = 1. Then pk > Pk+i for all k> 1. 

Proof Put 

OO 

g( S ) = {S- 1 )/(s) = -po + y(p k -l ~ Pk)s k 


g(s) — s — 1 + (1 — s) 2 [a + c(l — s) e ] l ^ 9 , 

g'(s) = 1 + c(l — s) e+1 [a + c(l — s) 9 ] -1 ^ 1 ^ 9 — 2(1 — s)[a + c(l — s) 6 ]~ lA 
= c(l-/( S )) 1+e + 2/( S )-l, 
g"(s) = (2-c(l + 8)(l-f(s)) e )ns), 
we see that g"(s) > 0, since 

G(s) := 2 - c(l + 0)(1 - f(s)) 9 > 2 - c(l + 0)(1 - Po ) 9 = 2- > 0 . 

a + c 

Furthermore, 

G , (s) = cd(l + e)(l-f(s)) 9 - 1 f , (s) 
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is absolutely monotone (as a product of two absolutely monotone functions), implying that g"{s) is absolutely 
monotone, so that 

k{k - l)(pk-i ~Pk) > 0, k > 2. 


4 Basic properties of f £ Ge 

In this section we distinguish among nine cases inside the collection of families {Ge}~i< 9 <i and summarize 
the following basic fomulas: f n (s), /(1), /'(1), /"(1). In all cases, except Case 1, we have a = f'(q). The 
following definition, cf [4], explains an intimate relationship between the Cases 3-5 with A = 1 and the Cases 
7-9 with A > 1. 

Definition 6 Let A > 1 and a probability generating function f be such that f(A) < A. We call 

/(») :=^ = f>.4*-V 

fc =0 

the dual generating function for f and denote q = qA _1 , so that f(q) = q. Clearly, f'(q) = /'((?). 

Case 1: 9 £ (0,1], a £ (1, oo), 

fn(s) = 1 - [a n (l - s)~ 9 + (a n - 1 )d]~ 1/e , d £ (0, oo). 

The corresponding theta-branching process is subcritical with m = a~ x ^ e . If 8 £ (0,1), then /"(1) = oo and 
for 9 = 1 we have /"(1) = 2 (a — l)a -2 d 

Case 2: 8 £ (0,1], a = 1, 

f n (s) = 1 - [(1 - s)~ 6 + nc\~ 1/e , c £ (0, oo). 

The corresponding theta-branching process is critical with either finite or infinite variance. If 9 £ (0,1), 
then /"(1) = oo and for 9 = 1 we have /"(1) = 2c. This is the only critical case in the whole family of 
theta-branching process. 

Case 3: 9 £ (0,1], a £ (0,1), 

fn(s) = 1 - [a"(l - s)~ e + (1 - 0(1 - g)- e ] ” 1/B , q £ [0,1). 

The corresponding theta-branching process is supercritical with m = a -1 / e . If 9 £ (0,1), then /"(1) = oo, 
and for 9 = 1 we have /"(1) = 2a _2 (l — a)(l — q)~ l ■ 

Case 4: 8 = 0, a £ (0,1), 

fn(s) = 1 - (1 - g) 1_Q ”(l - s) Q ”, q £ [0,1). 

The theta-branching process is regular supercritical with infinite mean. 

Case 5: 9 £ (—1, 0), a £ (0,1), 

fn(s) = 1 - [a"( 1 - S )' e l + (1 - a n )( 1 - g)l e l] 1/|fl| , q £ [0,1). 

The theta-branching process is non-regular with a positive 

l-/(l) = (l-a) 1 /l e l(l-g) 
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and infinite /'(1). 


Case 6: 9 = —1, a £ (0,1), 

fn(s) = a n s +(1-a n )q, ?e[0,1]. 

If q — 1, then the theta-branching process becomes a pure death process with mean m = a and /"(1) = 0. 
If q < 1, then the theta-branching process is non-regular with a positive 

1-/(1) = (l-a)(l ~q), 

/'(l) = « and /"(1) = 0. 

Case 7: d£ (0,1], a £ (0,1), A > 1, 

Ms) = A- [a n (A - s)~ e + (1 - a n )(A - q )- e ]- 1/e , q £ [0,1], 

If q = 1, then the corresponding tlreta-branching process is subcritical with the offspring mean m = a and 

f"(l) = (l+9)a(l-a)(A-l)- 1 . 

If q £ [0,1), the theta-branching process is non-regular with a positive 

1 - /(1) = (A - 1 )([a + (1 - a)(A - q)~ s (A - l) e ]~ 1/6 - 1), 

and 

/'(l) = o[a + (1 - a){A - g)- e (A - l)®]- 1 /*" 1 , 

/"(i) = a+0)0(1 - o)(^ - 9 rV - i) e ->+(i - o)(^ - 9 )- e (^ - i) e r 1/e ~ 2 - 

We have /(A) = A, and the dual generating function has the form of the Case 3: 

f(s) = 1 ~ [o(l - s)- e + (1 - o)(l - q)- 9 ]- 1 /*. 

Case 8: 6 = 0, a £ (0,1), A > 1, 

fn(s) = A — (A - q) 1 ~ an (A- s) a ”, q £ [0,1]. 

If q = 1, the theta-branching process is subcritical with the offspring mean m = a and 

/"(l) = a(l-a)(A-l)- 1 . 

If q £ [0,1), the theta-branching process is non-regular with a positive 

1 - f(l) = (A- qy- a (A - l) a -(A- 1), 

and 

f'(l) = ci(A-q) 1 - a (A-l) a - 1 , 
f"(l) = a(l-a)(A~q) 1 ~ a (A-l) a - 2 . 

We have f(A) = A , and the dual generating function belongs to the Case 4: 

f(s) = l-(l-q) 1 ~ a (l-s) a . 

Case 9: 8 £ (—1, 0), a £ (0,1), A > 1, 

M s ) = A — [a n (A-s)W + (l~a n )(A-q)W] im , «e[0,l]. 

If q = 1, then the theta-branching process is subcritical with the offspring mean m = a and 
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If q € [0,1), the theta-branching process is non-regular with a positive 

1 - /(1) = [a(A - l)l»l + (1 - a) (A - qf^ -(A- 1), 

and 

/'(l) = a[a + (1 - a) (A - q)W(A - l)-l*l] 1 /l*l- 1 e (0,1), 

/"(1) = (1 - |0|)o(l - a){A - qf\A - l)-!*!"^ + (1 - a)(A - q) m (A - s ) 

With 

f(A)=A-(l-a) 1 W(A-q)e(q,A), 

the dual generating function takes the form of the Case 5: 

/(s) = 1 - [a(l - s) |e| + (1 - a)(l - g) |e| ] 1/|f?l . 

5 Extinction and explosion times 

Recall that T = To A T\, and in the regular case T = Tq. In the non-regular case, when /(1) < 1, from 

P{n < T 0 < oo) = q - /„(0), 

P{n < Ti < oo) = /„( 1) - q, 

we obtain 

P(n<T< oo) = f n { 1) - /„(0). 

For our special family of branching processes we compute explicitly the distribution functions of the times 
T 0 ,7i,T. 

Cases 1-4. In these regular cases we are interested only in the extinction time: 

a-"/®[l + d-da- n ]-V®, Case 1, 

(1 + cn) -1 / 0 , Case 2, 

(1 - g)([l - a"(l - (1 - g)*)]- 1 / 0 - 1), Case 3, 

(i - g)[(l - g)- a " - i], Case 4. 

Cases 5, 7, 9. In these cases 

P(n < T 0 < oo) = (A - g)([l - a"(l - (A - ^A"*)]" 1 /® - 1), 

P(n < T x < oo) = (A - g)(l - [1 - a"(l - (A - g)®(A - I)-*)]" 1 /*), 

P(n < T < oo) = (A - g)| [1 - a"(l - (A - q) e A"*)]- 1 /® - [1 - a n (l - (A - g)®(A - l) _e )] _1/e }- 
Case 6. In this trivial case 

P(n < To < oo) = a n q , P(n < Ti < oo) = a n (l — g), P(n <T < oo) = a", 
and for q £ (0,1), 

1 


P(n <T 0 < oo) = < 


E(T 0 \T 0 < oo) = TOlTi < oo) = P(T) = 


1 — a 


Case 8. In this case 


P(n < T 0 < oo) = (A - g)[(A - g)- a "A a ” - 1], 

P(n < Ti < oo) = (A - g)[l - (A - q )~ an (A - 1)°"], 
P(n < T < oo) = (A - g) 1 - 0 " [A a " - (A - l) a "]. 
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Theorem 7 Consider a theta-branching process with 9 £ (—1,0] and A > 1. Let 9 —» 0 and A —> 1 in such 
a way that 

|0| ‘ log J3T V ’ re [ 0 ’°°]' 

Then for any fixed a £ (0, 1), q £ [0, 1), and y £ (— oo, oo), 

lim P(Tr - log Q e < y\Ti < oo) = e~ waV , 

e ->0 

where 

I |0|, r £ (0, oo], f 1, r £ {0} U {oo}, 

[ (log^zr) -1 , r = 0 > [ 1 - e_r > r-e(0,oo). 

The limit is a Gumbel distribution with mean lo f w ~~( where 7 is the Euler-Mascheroni constant. 

log a 7 ' 

Proof In view of 

P{Ti < n\Ti < 00 ) = [1 - a ra (l - {A - l)' e l(A - g)"l®l )] 1,m - 

it suffices to verify that 

[1 - ea y { 1 - (A - 1)1®!)] l,W -> e~ waV . 

Indeed, if r = 00 , then (A — 1)I* 9 ! -» 0, and 

[l - |0|a y (1 - (A - l) |e| )] 1/|S| ->• e~ aV . 

If r £ (0, 00 ), then {A — l)l e l —> e~ r , and 

[1 - \9\a y {l - (A - 1)1®!)] l,m -> e- ay(1 ~ e ~ r \ 

Finally, if r = 0, then 

l_(A-l)l 0 l~|0|/ e , 

and therefore 

[1 - ea y (l -(A- l ) 1 ® 1 )] 1/101 ->■ e~ aV . 


Corollary 8 If A = 1 and 9 £ (—1, 0), then for any fixed a £ (0,1) and q £ [0,1), 

lim P(Ti - log a |0| < y|Ti < 00 ) = e _a ”, y £ (- 00 , 00 ), 
6 —^0 

7/0 = 0 and A = 1 + e~ x P, e > 0, then for any fixed a £ (0,1) and q £ [0,1), 

lim P(Ti - log a e < y|Ti < 00 ) = e~ aV , y £ (- 00 , 00 ). 


6 The Q-process 

As explained in Ch 1.14, (I], for a regular Galton-Watson process with transition probabilities P n (i,j ), one 
can define another Markov chain with transition probabilities 


Qnifi j) 


jtf l P n (iij) 

7 n i 
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i> 1 , ,3 > 1 , 







where 7 = f'(q). The new chain is called the Q-process, and from 


i> 1 


7 n iq l ds 


( fn(sq )) 


fn(sq) _ / fnjsq ) y ~ 1 

fM ^ q > 


we see that the Q~P rocess is a Galton-Watson process with the dual reproduction and an eternal 

particle generating a random number n of ordinary particles with E(s K ) = 7 ^ 7 , see 0 . The Q-process in 
the regular case is interpreted in [T] as the original branching process ’’conditioned on not being extinct in 
the distant future and on being extinct in the even more distant future”. 

Exactly the same definition of the Q-process makes sense in the non-regular case, only now the last 
interpretation should be based on the absorption time T rather than on the extinction time To. Indeed, 
writing Pj(-) = P{-\Z 0 = j) we get for j > 1, 


and therefore, 

Pi(Zi =ji,...,Z n = j n \T > n+k) = Pi(Zi = j u ...,Z n = 

Jn+k\ L ) /n+fc! U l 

In the non-regular case, as k —> 00 we have /&(0) q and fk( 1) —> q. Thus, repeating the key argument of 

Ch 1.14, [T; for the derivation of the Q-process, 

Pi{Zi =ji,...,Z n = j n \T > n + k) -» P i (Z 1 = ji ,..., Z n = jn) j 7 

7 n iq l 

we arrive in the limit to a Markov chain with the transition probabilities Qn(hj)- 
By Theorem 3 from Ch. 1.11 in ill, 

7 ~ n Pn{i,j) iq^Vj, i,j > 1, 

where Q(s) = satisfies 

Q(f{s)) = lQ(s), Q(q) = 0 . 

In the critical case as well as in the subcritical case with Y^kL 2 Pkk log k = 00 the solution is trivial: Q(s) = 0. 
Otherwise, Q(s) is uniquely defined by the above equation with an extra condition Q'(q) = 1, so that the 
Q-process has a stationary distribution given by 

Qn(i,j) -t jq j l Vj, 

with 

^2jq J ~ 1 k > jS 3 = sQ'(sq). 

i> 1 

These facts concerning Q{s) remain valid even in the non-regular case. It is easy see from © that for 
our family with 0^0 and A > q, the generating function 

Q(s) = (A- s )~ e -(A~ q )~ e , 

is determined by parameters (0, A) and is independent of a = 7 . Similarly, for 0 = 0 we have 

Q(s) = log —r —-. 

A-q 

This leaves us with two cases when A = q = 1. In the critical Case 2 the answer is trivial: Q{s) = 0. In the 
subcritical Case 1, we have 7 = a -1 / 0 and 

(l-/( S ))- e + d = 7 - e [(l-s)- e + d], 
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which yields 


Q(s) = [(l-s)- e + d\- 1 / e . 

From these calculations it follows, in particular, that for our family of branching processes, in all subcritical 
cases, the classical a; log a; moment condition holds: 

OO 

^^Pkklogk < oo. 
fc =2 

Using these explicit formulas for Q(s) we can easily find the conditional probability distributions 


lim P(Z n 

n—too 


j\T >n)= bj, j > 1. 


For all cases, except the critical Case 2, we have 


b 3 s3 

3> 1 


Q(sq) 

Q( 0) ' 


Turning to the Case 2, recall that for any critical Galton-Watson process, there exists a limit probability 
distribution 


lim P(Z n = j\T 0 


n+ 1) = Wj, j > 1, 


such that 


Since 


J2 w 3 sJ 

3> i 


lim 


fn(sp 0 ) - /„( 0 ) 
fn(po) ~ fn( 0 ) ' 


Mspo) = 1 - [(1 - a(l - [1 + c]- 1 ' 9 ))- 0 + no}- 1 ' 9 , 

we obtain 

^ • [1 — s(l — [1 + c] -1 / e )]- e — 1 

2 ^ w 3 s = -"-■ 


7 Embedding into continuous time branching processes 

Recall that a Galton-Watson processes with generating functions /„ is called embeddable , if there is a semi¬ 
group of probability generating functions 

Ft+u(s) = F t (F u (s)), t £ [0, oo), u £ [0, oo), (3) 

such that / ra (s) = F n (s), n = 1,2,_ Although not every Galton-Watson process is embeddable, see Ch. 

III.6 in [Tj, in this section we demonstrate that all theta-branching processes are embeddable. 

Behind each semigroup ([3J there is a continuous time Markov branching process with particles having 
exponential life lengths with parameter, say, A. Each particle at the moment of death is replaced by a 
random number of new particles having a probability generating function 


h(s) = h 0 + h 2 s 2 + h 3 s 3 + - 


For such a continuous time branching process (Z t )ts[ o,oo) the probability generating function Ft(s) 
satisfies 

^ x 
h(x) — x 



Es Zt 


(4) 


(see [8j for a recent account of continuous time Markov branching processes). Our task for this section is for 
each / £ Qe to find a pair (h, A) such that f(s) = Ei(s). We will denote by p — ^hk the corresponding 

offspring mean number and by q the minimal nonnegative root of the equation h(s) = s which gives the 
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extinction probability of the continuous time branching process. 
Cases 1-3. For a pair 9 £ (0,1] and /i £ (0,1 + d -1 ], put 


h(s) = 1 - n(l - s) + 1 - s) 1+S - 


Taking successive derivatives of h it easy to see that it is a probability generating function with h'( 0) = 0. 
Next we show that using this h as the offspring probability generating function for the continuous time 
branching process we can recover f(s) for the theta-branching processes as -Fi(s) by choosing /i and A 
adapted to Cases 1- 3. 

Case 1. For a given pair a £ (0,1) and d £ (0, oo), put 


M = 


(1 + 8)d 


(l + 8)d+l’ 

In this subcritical case, applying |4]) we obtain for s £ [0,1) 


A = [(1 + 8 l )d + 9 i ]lna. 


rFt(s) 

A t = 

J S 

yielding the desired formula 


dx 


' ,Ft(s) d logj^ 


(1-0X1-*) +1^(1 -xy+t 


1-0+ T^e ei °s( 1 - x )’ 


, -i/fl 


F t (s) = 1 - (a*(l - s)~ e + (a 4 - l)d)' 

Case 2. For a given c £ (0, oo), put n = 1 and A = (1 + <? - 1 )c. Then by dH), we get 

F t (s) = 1 - ((1 - s)~ 9 + ctj 

Case 3. If /i > 1, then q = 1 — (— TMi±C)i /0 anc j ^} ie proposed h can be rewritten as 


h(s) = s + 


(l_ s) i+e_(l_ 9 )8(l_ s ) 


1 + 0- (l~q) e 

For a given pair a £ (0,1) and q £ [0,1) choosing 

A = [(1 + # -1 )(l — q)~ e — # -1 ] lna“ 


and applying d4j , we obtain 

F t (s) = 1 - [a\ 1 - s)- 9 + (1 - a*)(l - q)- 9 ]- 1 ' 9 . 

It is easy to see that f(s) = Fi(s) covers the whole subfamily Qg corresponding to the Cases 1-3. 

Notice that if 9 = 1, then h(s) = 1 — ^ ^ s 2 generates the linear birth and death process with h"( 1) = /.i. 

If 9 £ (0,1), then h"{ 1) = oo. 


Case 4. Consider a supercritical reproduction law with infinite mean 


h(s) 


s + (1 


ln(l - s) - ln(l - q) 

S > i _ l n( i _ q) 


For ho £ [0,1) this can be rewritten as 


OO fc 

!.(») = ho + (1 - Ao) EtfFTl)- 
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In this form with ho = 0, the generating function h appeared in [5] as the reproduction law of an immortal 
branching process. Earlier in [101 . this reproduction law was introduced as 


h(s) = 1 - (1 - h 0 )( 1 - s)(l - ln(l - s)). 

To see that the theta-branching process in the Case 4 is embeddable into the Markov branching process with 
the above mentioned reproduction law, use the first representation of h and apply Q. As a result we obtain 
for s ^ q, 

A t _ [ Ft(s) _ dx _ _ r Ft{s) ln(l - x) 

1 — ln(l — q) J s (l-z)(ln(l- 2 ;)-ln(l-g)) J s ln(l - q) - ln(l - a;) 

= ln[ln(l — s) — ln(l — g)] — ln[ln(l — F t (s)) — ln(l — g)]. 

Putting A = (1 — ln(l — g)) In a -1 , we derive 

Ft.{s) = 1 — (1 — g) 1_a *(l — s) a *. 


Cases 5, 7, 9. In these three cases the corresponding h and A are given by an extension of the formulas 
for the Case 3: 


h(s) = s + 


(A — s) 1+e — (A — q) e (A — s) 
(1 + 6) A s -(A- q) 6 


X = [(l + d- 1 )A e (A-q)- e 


0 _ 1 ]lna _1 . 


Turning to Definition [ 6 ] we see that this h in the Case 7 is dual to the h in the Case 3, and in the Case 9 it 
is dual to that of the Case 5. 


Case 6. In this trivial case the corresponding continuous time branching process is a simple death- 
explosion process with h(s) = q and A = In a -1 . 


Case 8. Similarly to the Case 4 we find that the pair 


h(s) 


s+ {A 


s) 


ln(A — s) — ln(A — q) 
l + lnA-ln(A-g) ’ 


A = (1 + In A — ln(A — g)) In a 1 , 


lead to 

F t (s) = A-(A-q) 1 - a \A-s) a ' 
Observe that this h is dual to that of the Case 4. 
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